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Abstract. For a dense and strongly interacting system, such as a nucleus or a strongly-coupled 
quark-gluon plasma, the foundation of hydrodynamics can be better found in the quantum descrip- 
tion of constituents moving in the strong mean fields generated by all other particles. Using the 
result that the Schrodinger equation and the Klein-Gordon equation can be written in hydrody- 
namical forms, we find that the probability currents of the many-body system in the mean-field 
description obey a hydrodynamical equation with stress tensors arising from many contributions: 
quantum effects, mean-field interactions, and thermal fluctuations. The influence of various contri- 
butions to the hydrodynamical motion is expected to vary with the temperature, as the quantum and 
mean-field stress tensors playing more important roles at low and moderate temperatures. 
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INTRODUCTION 

The foundation of hydrodynamics is usually presented within the kinetic theory, in 
which particles are described as a gas with weak interactions [1]. For a dense and 
strongly interacting system, such as a nucleus or a strongly-coupled quark-gluon plasma, 
the foundation can be better found in a quantum description of a system with constituents 
moving in the strong mean fields generated by all other particles [2, 3]. As the time- 
dependent Schrodinger equation can be cast into a hydrodynamical form Jii |2|], the 
evolution of the probability currents in the mean-field potential provides support for 
the validity of the liquid-drop model of a nucleus [3], as in Bohr and Wheeler f^, 
and the validity of the dynamical model of nuclear fission, as in Hill and Wheeler L^. 
The deviation of classical hydrodynamics from the quantum description is embodied in 

the presence of the quantum stress tensor pf) [0, [3|]. Quantum shell effects manifest 
themselves as nuclear shell effects superimposed on a smooth hydrodynamical liquid- 
drop background, and they lead to the intrinsic deformation in many nuclei [jTj]. 

In particle and nuclear collisions at high energies, a phenomenological relativistic hy- 
drodynamical description of the process is a reasonable concept, as pioneered by Landau 
dsl] and supported by experimental findings ^,[io, H ]. In the related search for the foun- 



dation of relativistic hydrodynamics using the Klein-Gordon equation, we encounter the 
difficulty that the naive probability density constructed from the wave function of the 
Klein-Gordon equation may be negative. The presence of a negative probability den- 
sity may appear to preclude its description in relativistic hydrodynamics. However, as 
pointed out by Feshbach and Villars [12], the negative probability density arises because 
the Klein-Gordon equation describes both particle and antiparticle degrees of freedom. 



Positive probability densities can be formulated when we separate out the particle and 
antiparticle degrees of freedom to write the Klein-Gordon equation as a set of two cou- 
pled time-dependent Schrodinger equations. We shall follow and modify the formalism 
of Feshbach and Villars [|l2fl in order to search for the foundation for relativistic hydro- 
dynamics in strongly-coupled systems. 



NON-RELATIVISTIC HYDRODYNAMICAL DESCRIPTION OF A 

NUCLEUS 

We begin by briefly summarizing the non-relativistic hydrodynamical description of a 
nucleus with strongly interacting constituents. As the Schrodinger equation for a single- 
particle in an external field can be cast into the form of a hydrodynamical equation 
a system of non-relativistic single particles interacting in their own mean fields can be 
described by non-relativistic hydrodynamics 10, [sj]- We can represent the single-particle 
state wave function for a state a in terms of the amplitude function 0^ and the phase 
function Sa as 



Mfairj) = ^airj)exp{-iSair,t)}. (1) 

In the time-dependent mean-field description we obtain the Euler equation for the 
probability current nu' 

where n = Y^a^a^I^ = La'^a^^'^iSa/'La'^a^I ^OT i = 1,2,3, is the occupation 
number for the state a. 
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-^pfj^ (r, t) = n(r, ^ y (Pr2n{r2,t)v{r, Y2) , (5) 

and v(r,r2) is the two-body interaction that generates the mean field. The stress tensor 
due to the mean-field interaction can also be given as 

^P\/M=n[r,)^{^^]^ (6) 



where W^^'^ is the energy per particle arising from the mean-field interaction. The mean- 
field stress tensor p^^^ is given explicitly by 

Thus, the total pressure arises from many sources: (i) quantum stress tensor from 

quantum effects, (ii) the thermal stress tensor pfj from the deviation of the individual 

velocity fields from the local mean velocities, and the (iii) mean-field stress tensor pfj^ 
from the mutual interaction between fluid elements. For example, for a nucleus in which 
the nucleons interact with the Skyrm interaction, we have lH] 

and 

p'if^l{to + lt3nydij, (9) 

where for the Skyrm I nucleon-nucleon interaction, to = —1057 MeV/fm^ is the two- 
body interaction strength, and to = +14463 MeV/fm^ is the three-body interaction 
strength. 

A nucleus is a strongly-coupled system. The quantum and mean-field stress tensors 
are the dominant component for the nuclear fluid at low and moderate temperatures. 

The thermal stress tensor p^'j can take on different values, depending on the occupation 
numbers of the single-particle states that determines the degree of thermal equilibrium 
of the system. For a thermally equilibrated system, the thermal stress tensor is 
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which is small for low and moderate temperatures. We observe that when + | ^ 

pfj for a strongly interacting system at low and moderate temperatures, there can be 
situations in which the system can behave quasi-hydrodynamically even though the state 
of the system has not yet reach thermal equilibrium, as is evidenced by the presence of 
nuclear collective vibrational and rotational states at low and moderate temperatures. 
In this case, the hydrodynamical state is maintained essentially by the quantum stress 
tensor and the strong mean-field stress tensor, and not by the thermal stress tensor. 



RELATIVISTIC HYDRODYNAMICS AND THE KELIN-GORDON 

EQUATION 



In the collision of high-energy particles and nuclei, the relativistic hydrodynamics de- 
scription was pioneered by Landau Subsequent considerations for the evolution of 
quark-gluon matter in the work of Bjorken [l3l Baym et al. [14], OUitraut [SQ, and 
many others [.171 lead to successful investigations of the hydrodynamics of strongl y in - 
teracting matter in extreme conditions in relativistic heavy-ion collisions rt9L lld.lllLll7|l . 

In the search for the foundation of relativistic hydrodynamics using the Klein-Gordon 
equation, we can overcome the apparent difficulty of a possible negative naive probabil- 
ity density by re-writing the Klein-Gordon equation as a set of coupled time-dependent 
Schrodinger equations for particles and antiparticles, as in Feshbach and Villars [12]. 

A general solution of the Klein-Gordon equation contains two components that can 
be represented by a column vector, 

where x+ and x have positive norms, \x±?'^ which can be interpreted as the proba- 
bility densities of particles and antiparticles respectively. We find it necessary, however, 
to modify the formulation of Feshbach and Villars so as to obtain a set of coupled 
equations that facilitate the re- writing of the Klein-Gordon equation in the hydrodynam- 
ical form. We re- write the Klein-Gordon equation for the wave functions and X- as 



m 



[ihdt - e±A«)z±= 2^^^^ {(tV - e±A)2+ {m + S)^+ [{E - e,A^)^~ ihdt{E - e„AO)]| x± 



where E is the positive expectation value of the operator ihdt for the single-particle state, 
e± = ±e, Ca = Vflg, and is the conserved particle number of the state. 

KLEIN-GORDON EQUATION IN HYDRODYNAMICAL FORM 



To obtain the equation of continuity and the Euler equation from Eq. ([121) . we write the 
wave functions x± in terms of the amplitude and phase functions as in Eq. ([T]), 

X±{r,t) = (^)±(r,0exp{/5±(r,0}. (13) 
We construct X± x ([I2l)-Z± x dHj)* and obtain 

dt[{E - eaA'')(l)l] + V ■ [(j)liVS± - e±A)] = X±, (14) 



=F' 

(12) 



where 
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+ [{m + yY + {E- eaAf] {xlx^ - X±X^) 

+ [zM(£-e«A«)]Uiz;+Z±ZT)- 



(15) 



As the quantities X+ and X are not separately a full divergence, the total number of 
particles and antiparticles are not separately conserved. However, the difference of the 
particle number and antiparticle number of the two components satisfies the equation 

dt[{E-eaA''){^l-^^_)]+V ■[^l{VS+-e+A)]-V .[f_{VS--e-A)]=X+-X-.{\6) 

It can be shown that the difference (Z+ —X-) is a complete divergence. Therefore, the 
quantity 
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is a conserved quantity. 

Hydrodynamical description will be appropriate after the active pair production stage 
has passed and the expansion of the system is driven by a slowly varying mean field. 
It is this type of motion, with the suppression of pair production, for which we wish to 
provide a hydrodynamical description. Under such a circumstance, one can speak of a 
single-particle system with a definite particle number which can take on the value 
Vfl = 1 for a particle and Va = — 1 for an antiparticle. From Eq. ([T2)) . the equation for the 
phase function S± for this simplified case is 



(<9,5± + e±A0) = 
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(18) 



For this case with suppressed pair production, = Vae = e±. We take the gradient of 
the above for / = 1 , 2, 3, and multiply by 0^ [E — e±A^) . Using the equation of continuity 
for this simplifying case with the suppression of pair production, we obtain 



{E~e±A^)^l{ViS±-e±A') 
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We can identity the fluid energy density £± as 



(19) 
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as it corresponds to the fluid energy density for a fluid element at rest. The fluid is 
characterized by a relativistic 4-velocity u^. We can identify 



E - e±A 







and u 
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(21) 



which obeys (m±)^ — {u±Y = 1 in the absence of the quantum effects. We can write the 
equation of motion for the probability currents in terms of the hydrodynamical equation 



m 



7=1 



,7=1 



2(£-e±A0)(m + ^) 
-(£-e±A0)<^2^V,5±-e±A' 



(V2(^)±) /(^)± - ( V5± - ^^±A)2] - (m + ^)2 + (£ - e±A0)2l 



(m + ^)2 
£ \<l>l{'^jS± - e±A]){ViS± - e±A' 



m+y 



(m + ^)^ 



(22) 



where i = 1,2,3. This is the Klein-Gordon equation in hydrodynamical form for its 
particle and antiparticle probability currents in external fields for the case with the 
suppression of pair production. The first two terms on the left-hand side correspond 



to dp^T^ with = £±u^u'^ for /i = 0, 1,2,3. The third term on the left-hand side 
is the quantum stress tensor arising from the spatial variation of the amplitude of the 
single-particle wave function [|2D, 



P±ij 
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which comes from the quantum nature of the fluid. The first two terms on the right- 
hand side contain forces coming from the scalar interaction, the electric field F^^ and 
the magnetic field F'-', as in classical considerations. The third term on the right-hand 
side is the relativistic correction to the time-like part of the vector interaction, and the 



last two terms represent relativistic corrections associated with the spatial and temporal 
variation of the scalar interaction. Thus the dynamics of the probability currents obey a 
hydrodynamical equation, with forces on the fluid element arising from what one expects 
in physical and classical considerations. The additional element is the presence of the 

quantum stress tensor p^f that comes from the quantum nature of the fluid. 



RELATIVISTIC HYDRODYNAMICS FOR A MANY-BODY 
SYSTEM IN THE MEAN-FIELD DESCRIPTION 

A many-body system in the time-dependent mean-field description consists of a col- 
lection of independent particles moving in a time-dependent self-consistent mean-field 
potential generated by all other particles [0,[2L[3]- Each single-particle state \]/a is char- 
acterized by a state label a, particle number v, energy Ea, and occupation number Hq. 
For simplicity, we consider the case in which the mean-field potential arises from a 
scalar two-body interaction v.5(ri,r2) and a time-like vector interaction vo(ri,r2). We 
further neglect the last three terms on the right-hand side of Eq. (|22l) which represent 
higher-order relativistic corrections. The equation of motion for the energy density Sav 
and velocity fields for i = 1,2,3 in the single particle state a and particle type v is 
then 



3 3 

- -*ivV,-^ + ^^*^-^. (24) 

where, in the frame with fluid element at rest, 

^(r,0 = I dhj n(r2,0v.(r,r2) (25) 

and 

A°(r,0 = y^^r |n+(r2,Oe++n-(r2,Oe-|vo(r,r2). (26) 

Here n+ = E+"a+^a+' = E-"a-0a-' and n = We consider a strongly 

interacting system in which the number of particles and antiparticles are equal so that 
'^+{'^2) = "-(^i) and n+(r2)e+ + n_(r2)e_ is zero. Then the contribution from the 
second term on the right-hand side of Eq. (|24|) is zero. Multiply Eq. (|24|) by and 
sum over {a, v}, we get 

3 

av j=l av 

+^ i w!?,),,) + il,nav(^l)V,y = 0. (27) 



We define the total energy density e by 

Y.'^avSaV = e, (28) 



av 



and the average 4-veIocity u by 

M = J^^avCavWav/e- (29) 



av 



We can introduce the thermal stress tensor pfj' for {«, j} = 1,2,3 as the correlation of 
the deviations of the single-particle velocity fields from the average 

J^«aveav(wLv - "')("^v - ^■') = pfj- (30) 
av 

For the case with the suppression of pair production, we obtained the equation of 
hydrodynamics 



d,i8u%') + t { V, {eu^u^+p^ +p^) + ^yjpl]^ } = 0, (31) 



3 

I 

where the total quantum stress tensor is 



2 2 
Pif = -T-^'L«av0av5o + -L«avV,-0avV,-0av, (32) 



and the pressure due to the interaction pfj^ is 



^p\f{r,t) = n{r,t)Viy{r,t) = n{r,t)^ J dh2n{r2,t)vs{r,r2). (33) 



The mean-field stress tensor pjj^ can also be given as 



on 
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(34) 



where W^^^ is the energy per particle arising from the mean-field interaction. As an 
illustrative example, we can consider a density-dependent two-body interaction 

v,(r,r2) = [a2 + a3n((r + r2)/2)]5(r-r2). (35) 
The contribution of the mean-field interaction to the stress tensor is then 

pS;^ = ^[a2 + 2a3n(r)]n%-, (36) 



whose magnitude increases with the density and the strengths of the interaction. 

The thermal stress tensor p^lj can take on different values, depending on the occupa- 
tion numbers of the single-particle states that determine the degree of thermal equilib- 
rium of the system. In the time-dependent mean-field description, the motion of each 
particle state can be individually followed [|l9n . The occupation numbers of the 
single-particle states will remain unchanged, if there are no additional interaction be- 
tween the single particles due to the residual interactions. When particle residual inter- 
actions are allowed as in the extended time-dependent mean-field approximation [j20(], 
the occupation numbers will change and will approach an equilibrium distribution as 
time proceeds. 

We note that as in the non-relativistic case, the total pressure arises from many 
sources. We come again to the observation that in situations when \p)f + pf'}\ ^ 

pfj^ for a strongly-coupled system at low and moderate temperatures, there can be 
situations when the system behaves quasi-hydrodynamically, even though the state of 
the system has not yet reached thermal equilibrium. In this case, the hydrodynamical 
state is maintained mainly by the quantum stress tensor and the strong mean fields. 



SUMMARY AND DISCUSSIONS 

For dense systems with strongly interacting constituents, a reasonable description of 
the system can be formulated in terms of constituents moving in the strong mean fields 
generated by all other particles. From such an analysis, we find that the probability 
currents of the system obey a hydrodynamical equation with the stress tensor arising 
from many contributions. There is the quantum stress tensor that arises from quantum 
effects, there is the thermal stress tensor that arises from the deviation of the single- 
particle velocity fields from the average velocity fields, and there is the mean-field stress 
tensor that arises from the mean-field interactions. 

The importance of the three different contributions depend on the physical situations 
that are present in the system. For low temperature dynamics for which the quantum 
effects and mean fields are important, the dynamics of the strongly-coupled system will 
be dominated by the quantum stress tensor and the mean-field stress tensor. In that case, 
the degree of thermalization is of less significance as the thermal stress tensor plays a 
minor role. On the other hand, for very high temperatures for which the magnitude of 
the temperature far exceeds the strengths of the mean-field interactions and the quantum 
pressure, the thermal stress tensor plays the dominant role. The magnitude of the thermal 
stress tensor will depend sensitively on the degree of thermalization of the system. In 
between these limits, one can envisage the transition from the quantum and mean field 
dominating strongly-coupled regime to the thermal pressure dominating weakly-coupled 
regime as temperature increases. 

What we have discussed is only a theoretical framework that exhibits clearly the 
different sources of stress tensors. To study specifically the dynamics of the quark-gluon 
plasma for example, it will be necessary to investigation the specific nature of different 
constituents and their interactions. Nevertheless, the general roles played by the different 
components of stress tensors can still be a useful reminder on the importance of the 



quantum and mean-field stress tensors in the strongly-coupled regime, at temperature 
just above the transition temperature Tc. The dominance of the quantum and mean- 
field stress tensors over the thermal stress tensor may imply that hydrodynamics may 
commence at an early stage even when thermal equilibrium has not been reached, as the 
total stress tensor may be hardly affected by the variation of the thermal stress tensor. It 
will therefore be of great interest to study whether the fast onset of the hydrodynamical 
motion as suggested by the elliptic flow measurements [[l6l] may be caused by the 
dominance of the mean-field and quantum stress tensors for a strongly-coupled quark- 
gluon plasma. For T slightly greater than such a possibility may perhaps be the case. 
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